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ABSTRACT 

The demonstration of repeated gamma-ray bursts from an individual source would severely 
constrain burst source models. Recent reports (Quashnock & Lamb 1993; Wang & 
Lingenfelter 1993) of evidence for repetition in the first BATSE burst catalog have 
generated renewed interest in this issue. Here, we analyze the angular distribution of 
585 bursts of the second BATSE catalog (Meegan et al. 1994). We search for evidence 
of burst recurrence using the nearest and farthest neighbor statistic and the two-point 
angular correlation function. We find the data to be consistent with the hypothesis that 
burst sources do not repeat; however, a repeater fraction of up to about 20% of the observed 
bursts cannot be excluded. 

Subject headings: gamma-rays: bursts 
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1. INTRODUCTION 

The observed isotropic and inhomogeneous spatial distribution derived from BATSE 
burst data (Meegan et al. 1992; Fishman et al. 1994) severely constrains possible Galactic 
distributions and argues in favor of sources at cosmological distances. Although neutron 
stars in an extended Galactic halo were considered as an alternative to cosmological models, 
the observational constraints are now so severe that a halo origin of bursts appears unlikely 
(Hakkila et al. 1994; Hartmann et al. 1994b; Briggs et al. 1995). In contrast, cosmological 
models naturally explain the observed isotropy and inhomogeneity, but some of these are 
incapable of producing multiple bursts from one source. 

The absence of an excess of overlapping error circles in pre-BATSE burst localizations 
provides a model-dependent upper limit of ~10 years on the burst repetition time scale 
(Schaefer & Cline 1985; Atteia et al. 1987). Several recent reports, however, have cited 
evidence for repetition on much shorter times. Quashnock & Lamb find an excess of close 
neighbors in the IB burst catalog, and conclude that a large fraction of classical bursts 
repeat on timescales of order months (Quashnock & Lamb 1993). Wang & Lingenfelter 
(1993) claim evidence for repetition with recurrence times perhaps as short as days from 
one particular location (0855—00). Both Quashnock & Lamb and Wang & Lingenfelter 
interpret their results as evidence for multiple repetitions, of about five observed bursts per 
repeating source. Also, the coincidence of the locations of GRB940301 and GRB930704 
determined with COMPTEL has a 3% chance probability, suggesting a possible repeater 
(Kippen et al. 1995). 

To investigate burst repetition, we test the null hypothesis that "the angular 
distribution of GRBs is consistent with the isotropic distribution", i.e., equal probability 
per unit solid angle. Tests of isotropy have varying sensitivities to clustering, which can 
indicate the presence of repetition, as well as to large-scale anisotropics. In this paper we 
focus on the implications for burst repetition of the two-point correlation function and the 
nearest neighbor test. We do not consider possible time dependent repetition, as suggested 
by Wang & Lingenfelter (1993), which is discussed in a separate paper (Brainerd et al. 
1995). Other clustering tests are considered in Hartmann et al. (1995). We consider 
only the "classical" gamma-ray bursts, which are distinct from the Soft Gamma Repeaters 
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(Kouveliotou 1994). 

We analyze the 2B catalog of bursts observed by BATSE between 1991 April 19 and 
1993 March 9, comprising 585 bursts (Meegan et al. 1994) and various subsets. Data after 
1992 March contain numerous gaps due to CGRO tape recorder errors. For 100 bursts 
that were most seriously affected by these gaps, the locations are determined with MAXBC 
data, which consists of the background-subtracted maximum rates in each detector on a 
1.024 second timescale in the 50 to 300 keV energy range. Koshut et al. (1994) find that 
the total location error (systematic + statistical) is about 7° for bright bursts located 
using MAXBC data. Since this is larger than the usual systematic error of about 4°, the 
inclusion of MAXBC-located bursts might obscure a repeater signal. 

Our subsets, which are listed in Table 1, are of consecutive triggers limited by the 
listed trigger numbers. We apply two cuts to some of the data sets: data sets marked with 
"yes" in the MAXBC column contain bursts located with the MAXBC data, while data 
sets with "no" in this column have such bursts removed. Similarly, data sets with "yes" 
in the column "> 9°" contain bursts with statistical location errors > 9°, while data sets 
with "no" in this column do not. The total location error of a burst is estimated as the 
rms sum of its statistical error and a 4° systematic error. 

We define / as the fraction of all observed bursts that can be labeled as repeaters and 
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v as the average number of observed events per source observed to repeat (thus v > 2). 
We define Nb as the number of observed bursts and n r as the number of sources from 
which two or more bursts were observed. These quantities are related via n r = fNs/v- 

2. TWO-POINT ANGULAR CORRELATION FUNCTION 

One mathematical function that is used to test for anisotropy is the two-point angular 
correlation function, w(8) } defined in the following manner: for an ensemble of points 
distributed on the sky, the average number of pairs with angular separation 6 within the 
solid angle dQ, is (e.g., Peebles 1980) 

Nr -1 

dN p = — [l + w(6)]dn. (1) 

The application of angular correlation analysis to GRB data was introduced by Hartmann 
& Blumenthal (1989) for point sources and refined for fuzzy sources by Hartmann, Linder, 
& Blumenthal (1991). If some of the observed bursts are repeaters, each of their positions 
on the sky will be displaced from the source position by a distance of order the total 
location error. If the location error is characterized by a Gaussian distribution with a 
standard deviation of cr* <C 1, then the the observed correlation function is 



(N B -l)w(0) = f(u-l) 



2 / e , 



(2) 



This equation shows that excess correlation is spread over an angular scale of cr* (~ 7° for 
BATSE). The negative correlation at larger angles occurs because the correlation function 
must integrate to zero. For a given repeater fraction / the strength of a repetition signal 
in the data increases with v — 1. 

Figure 1 shows the two-point correlation functions for data sets 1, 5 and 6 (Table 
1), which are the revised IB catalog, the 2B catalog and the 2B catalog less MAXBC- 
located bursts. In the IB data two regions with excess at the 2a level are apparent, 
one near 0°, which is interpreted in Quashnock & Lamb (1993) as evidence of repetition, 
and one near 180°, which Narayan & Piran (1993) use to argue against the repetition 
interpretation (but see Quashnock & Lamb 1994). Here we use the final localizations of 
the 2B BATSE catalog and find that there are no significant angular correlations on any 
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Fig. 1 — The angular correlation function of gamma-ray bursts. Shown are the 
results for 262 bursts in the IB catalog (data set 1), the full 2B set of 585 bursts 
(data set 5), and the modified 2B set with 485 bursts in which MAXBC events 
(see text) were removed from the sample (data set 6). The addition of second 
year data clearly reduces both excesses near 0° and near 180° originally found in 
the IB set. 
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scale in the current data set. Burst data obtained after the IB period do not show the 
excesses that are apparent in the IB set. The combined data sets still show some residual 
effects of the IB excesses but the deviations are not significant. Consistency with the null 
hypothesis of zero correlations is evaluated with the Kuiper (1960) statistic, which has 
certain advantages over the usual Kolmogorov-Smirnov (KS) test in the current context 
(Hartmann et al. 1995). We have also applied various cuts in angular resolution, using 
subsets of the data with better localization accuracy. None of these sets shows evidence 
for significant deviations from zero clustering. To derive upper limits on the presence of 
observed repeaters, we fit the data with the model correlation function (eq. 2). We set 
cr* = 7°, the mean positional uncertainty of BATSE locations, and fit the first two 4° 
wide bins of the correlation function using various values of f(v — 1). The fits become 
unacceptable at the 99% confidence level if f(v — 1) exceeds 0.2. The most difficult 
case to detect is v = 2: the limit for this case is / ^ 0.2. Repetition has been previously 
reported for v ~ 5: the limit for this case is / ^ 0.05. 



Another approach to detecting burst recurrences is the nearest neighbor test, which 
tests whether the separations between bursts are consistent with the separations found for 
the isotropic distribution. For isotropically distributed bursts one expects the cumulative 
distribution of nearest neighbors to be (Scott & Tout 1989) 



Burst repetition will create small-scale anisotropics in the burst densities. The nearest 
neighbor test can indicate the existence of such anisotropics if the average distance between 
bursts is greater than the location error. For BATSE, this requires that the sample size 
be less than about 500 bursts (Brainerd et al. 1994), although multiple repetitions may 
be detected in larger samples. 

The First BATSE Gamma-Ray Burst Catalog was analyzed by Quashnock & Lamb 
(1993) for burst repetition by comparing through the KS statistic the cumulative 
distribution of nearest neighbor separations with that expected for a uniform sky 



3. NEAREST NEIGHBOR TEST 
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Fig. 2 — Nearest neighbor cumulative distributions from the 2B catalog, plotted 
as functions of 1 — cos 6, where 6 is the angle to the nearest neighbor. The data are 
plotted as a histogram while the model curve for isotropy is plotted as a smooth 
curve. The four plots are for data sets 1, 3, 7, and 8 of Table 1. Data sets 1 and 7 
correspond to the IB catalog while data sets 3 and 8 correspond to the 2B — IB 
catalog. Data sets 1 and 3 are consecutive sets of 262 gamma-ray bursts. Data 
sets 7 and 8 are bursts with position errors < 9°. 



distribution. This was done for the full catalog of 260 bursts and for various subsets. They 
found a deviation from isotropy of 2% significance for the full catalog and of 1.1 X 10 -4 
significance for the 202 bursts with statistical errors less than 9 degrees. The selection of 
the 9 degree error cut maximizes the signal but introduces uncertainty in the calculation 
of statistical significance, since the value of 9 degrees was not specified a priori. Such 
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techniques are useful for exploring a data set for unanticipated effects but must be treated 
as predictions for subsequent data sets. 

Results of our nearest neighbor analysis of the Second BATSE Catalog and various 
subsets are given in Table 1. Figure 2 shows the nearest neighbor cumulative distribution 
for four of these subsets. We find the maximum deviation D of each data set from the 
isotropic cumulative distribution and derive the KS statistic K = D^JNb-, where Nb is 
the sample size. The significance S of the magnitude of K — that is, the fraction of trials 
that produce a greater deviation from the model curve — is determined through Monte 
Carlo simulation; the usual analytic formula is invalid for this analysis because the nearest 
neighbors are not statistically independent. The results are given in Table 1 for both the 
celestial coordinate frame (A' ce i and S ce \) and the CGRO coordinate frame (K gro and S gro ). 
The analysis in CGRO coordinates is particularly sensitive to systematic effects relating 
to the angular response of the BATSE detectors. Such effects would be less apparent in 
the celestial coordinate frame because the CGRO orientation is routinely changed at one 
or two week intervals. The effect seen by Quashnock & Lamb is reproduced in our data set 
7, while the remaining subsets exhibit no statistically significant deviation from isotropy. 

An upper limit on the number of repeating sources can be found from both the nearest 
neighbor test and the farthest neighbor test. Through Monte Carlo simulation we derived 
these limits for an isotropic distribution of burst sources. The model that the various data 
sets were tested against consists of n\ sources that each produce one observed burst and n r 
sources that each produce v observed bursts. The burst locations of the repeating sources 
are given a Gaussian distribution with a 9° standard deviation about the source location. 
Limits on the burst repeater fraction are derived from the maximum deviation of the 
various data sets from the average nearest and farthest neighbor cumulative distributions 
of the Monte Carlo simulations. Because of the computational demands of this analysis, 
we determine the significance of a deviation from the probability of a similar deviation in 
the no-repeater model. We have verified that this approximation is adequate by deriving 
some significances using simulations of repeater models. The upper limit on the repeater 
fraction for all sets is roughly the average value of the upper limit found for each data 
set, which is the smaller of the limits derived from the nearest and the farthest neighbor 



8 



statistic. We find from the Monte Carlo simulations that the nearest and farthest neighbor 
statistics limit the quantity f(v — l) 1 ' 2 jv. Except for data set 7, we find the limit to be 
f(u - iy- 2 /u < 0.2. For v = 2 the limit is / < 0.4, while for v = 5 it is / < 0.2. 

We test data set 8, which was selected with the 9° criterion of Quashnock & Lamb, 
against the null hypothesis and find no evidence of burst repetition. As explained below, 
we calculate that BATSE remains sensitive to repeaters, so the results of data sets 7 and 8 
are contradictory. Considering the difficulty of evaluating the significance of the evidence 
for repeaters found by retrospective analysis of the IB catalog and the non-confirmation 
of the effect in the post-IB data, we conclude that the data are consistent with the null 
hypothesis of isotropy. 

4. DISCUSSION 

If the MAXBC-located bursts are removed from the 2B catalog, the effective exposure 
decreases. Exposure is defined here as the fraction of bursts above trigger threshold that 
would be observed, and is less than unity due to earth blockage, SAA passes, and intervals 
during which the burst trigger is disabled. Our investigation of the effect of exposure on 
the detectability of repeaters finds that if the number of detectable bursts from a repeating 
source is large, then the average number of bursts observed for each source remains large, 
and the fraction of bursts identified as bursts from repeater sources changes little, but if 
the number of bursts from a repeater source is small, so that the average observed number 
of bursts is small, then the number of repeater sources that appear as single burst sources 
increases significantly as the efficiency of detecting a burst drops. The effect of this is to 
decrease the fraction of bursts that are identifiable as bursts from repeaters. 

Figure 3 illustrates the effect of a varying exposure on the detectability of repetitions 
in one specific model — an ensemble of sources that each produce 10 bursts above threshold, 
not all of which are detected. The solid line, referenced to the left axis, is the fraction 
of events that appear to have at least one companion burst. Our simulations have shown 
that the strength of the signal in the nearest neighbor test is approximately proportional 
to this curve. The dashed line, referenced to the right axis, is the average number of 
bursts observed from sources that produce at least two observed bursts. The strength 
of the signal in the two point correlation function is approximately proportional to this 
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Exposure 



Fig. 3 — Efficiency of observing burst repetition as a function of sky exposure, for 
a model in which each source produces 10 outbursts. The solid line, referenced 
to the left axis, is the fraction of events that can be identified as repeaters. The 
dashed line is the average number of bursts observed from sources that produce 
at least two observed bursts. The right hand vertical line is the exposure for 
the IB catalog and the left hand vertical line is the exposure for the 2B — IB 
catalog. 

curve. The vertical dashed line at 0.34 indicates the exposure of the IB catalog. Here, 
an average of about 3.5 bursts will be observed from each repeating source, and about 3% 
of the observed bursts will be misidentified as non-repeaters. Note that this number of 
observed repetitions per observed repeater is slightly less than the value v <; 4 suggested 
by Quashnock & Lamb (1993), based on the angular scale of the dumpings seen in the IB 
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catalog. A smaller source repetition rate than our choice of 10 produces a smaller observed 
repetition rate and a larger fraction of misidentifications. The vertical dashed line at 0.25 
indicates the exposure of the post-IB portion of the 2B catalog when MAXBC-located 
bursts have been removed. Here, an average of about 3.0 bursts will be observed from 
each repeating source, and about 8% of the observed bursts will be misidentified as non- 
repeaters. For this specific model, a change in exposure of only 26% has a negligible effect 
on the burst repetition limit derived for the 2B catalog from the nearest neighbor analysis. 
The limit from the two point angular correlation function is increased by ~ 15%. 

While we find the data to be consistent with no burst repetition, repetition at some 
level cannot be excluded. Using a simple model, we place upper-limits on repetition using 
the two-point angular correlation function and the nearest and farthest neighbor statistics. 
In all cases the two-point angular correlation function places a tighter limit on the fraction / 
of observed bursts that could be from repeating sources than the nearest neighbor statistic. 
The limit for v = 2 observed bursts per repeating source is / ^ 0.2 and for v = 5 the limit 
is / ^ 0.05. We do not place a limit on the fraction of sources that emit multiple bursts 
because such limits are highly model dependent. 

Several factors will improve these results. Statistical limits will be reduced as BATSE 
continues to accumulate burst locations. Flight software changes since the 2B catalog have 
eliminated the need for MAXBC locations. The current daily exposure exceeds that of 
the IB era, primarily because of reduced solar activity. Finally, we continue to refine the 
burst location algorithm to reduce systematic errors. As a consequence, analysis of the 
forthcoming 3B catalog will improve the constraints on burst repetition. 

Acknowledgements: This work is supported in part by NASA grant NAG 5-1578. 
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